One dimensional gapless magnons in a single anisotropic ferromagnetic nanolayer 
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Gapless magnons in a plane ferromagnet with normal axis anisotropy are shown to exist besides the 
usual gapped modes that affect spin dependent transport properties only above a finite temperature. 
These magnons are one dimensional objects, in the sense that they are localized inside the domain 
walls that form in the film. They may play an essential role in the spin dependent scattering 
processes even down to very low temperatures. 
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While charge and its transport play an essential role 
in carrying and processing information, spin is the origin 
of magnetism and is the basis to most of the storage de- 
vices in information-related technologies. There is much 
interest in combining charge and spin properties in or- 
der to expand the possibilities for applications into the 
so called spintronics 0. This field, that has started with 
the discovery of giant magnetoresistance (GMR) in al- 
ternating Fe/Cr layers Q, has rapidly become on focus 
with the report of magnetic materials with high transi- 
tion temperatures (Tc) such as Gai-^Mn^As (Tc ~ 110 
K) and Tii_ K Coz0 2 ( T c ~ 300 K) for instance. 
With the upscaling of Tc toward room temperatures re- 
alizations of practical ideas ensue, pointing in the direc- 
tion of spin-based devices analogous to p-n transistors 
II- Functioning principles of such devices rely on the 
existence of domain walls (DW) in finite ferromagnetic 
samples. The role of DW on transport of electrical cur- 
rent has been shown in the early experiments on iron 
whyskers^ and since then subjected to intense experi- 
mental and theoretical investigation Q. The possibility 
that GMR originates directly from DW may yield the 
fabrication of GMR devices from single layered mate- 
rials taking simultaneous advantage of spin and charge 
degrees of freedom. With the improvement of growing 
techniques, the possibilities of producing 2D ferromag- 
netic samples with DW of few tenths of A has become 
reality 0- 

In this Letter we report on a magnetic excitation that 
has not yet been accounted for the scattering of electrons 
by DW. One-dimensional (ID) gapless magnons prop- 
agating inside a DW of a single-layer ferromagnet are 
shown to exist in addition to the usually known gapped 
magnons. Starting from an appropriate model hamilto- 
nian, we solve it for the continuum limit in 2D and find 
that these modes are confined in the ID DW, as illus- 
trated in FigQ] Due to the absence of gap, the inelastic 
coupling of these modes to electrons should manifest at 
very low temperatures, as opposed to the usual magnon 
scattering processes that become important only at tem- 
peratures comparable to the anisotropy gap. We provide 
a theory to include the effect of the gapless magnons on 
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FIG. 1: Domain wall in 2D showing the ID channel at its 
center where gapless magnons are expected to propagate in 
the y direction; the amplitude of such a mode is depicted by 
the bell-like curve shown on top. 



conduction electrons. We then present results for the re- 
sistivity in a 2D system as a function of the DW size, 
the electronic density, and the temperature for realistic 
values of the parameters. 

In order to model a system of independent electrons 
interacting with local 1/2 spins, we will use the effec- 
tive hamiltonian H = JC + Hh + Hi, where JC stands for 
the electron's kinetic energy operator — h 2 W 2 /2m, while 
Tin accounts for the interaction between local spins (Si), 
described by a 2D ferromagnetic Heisenberg spin hamil- 
tonian 
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on a square lattice of constant a. Here A > gives 
the strength of the on-site anisotropy and J > is the 
exchange amplitude. Interaction between electrons and 
the local spins is considered through 



Hi = -J K a 2 Y,°- ■ SiS(n - r) 



(2) 



where Jk is the exchange integral and <j is the spin of 
the electron. 

First, we consider the static DW solution of Hh for a 
7r-domain. In the continuum the following hydrodynamic 
equations are obtained for the magnetization, 
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with j = JS 2 a 2 /h and A = a/V2A. The fields and 
parametrize the magnetization through the usual relation 
S(r) = S (sin 9 cos (j), sin 9 sin (j), cos 9). Such a hydrody- 
namic treatment suits the low energy excitations studied 
in what follows, hence the continuum limit is not only 
appropriate but becomes exact for the long wavelength 
phenomena handled below. 

It is possible to find a static solution for these coupled 
equations in the form of a it DW. Assuming <f> = <j) is con- 
stant, 9 = 9(x) (see Fig^l, and using the boundary con- 
ditions linij.^zpoo cos 9(x) — ±1, the coupled equations 
can be easily integrated to give S(r) = S[x sech(a;/A) — 
ztanh(x/A)] if we conveniently choose (f> = 0. 

Hence a ir DW like the one illustrated in Figgis admit- 
ted as an equilibrium configuration in the 2D ferromag- 
netic system described by Tin- One sees that A measures 
the size of the spin-inversion region thus giving the size 
of the DW. From the definition of A it is clear that the 
DW extends throughout the sample unless the anisotropy 
is different from zero. Strictly, the DW configuration is 
not the absolute minimum of energy in our model, al- 
though it is a stable solution. If one however takes the 
effect of dipolar interactions, as in a finite sample, the 
ground state becomes, in general, a DW distortion. Yet, 
if the temperature is lowered, and as long as the on-site 
anisotropy A > 0, one expects to cross a transition tem- 
perature below which there are no DW structures in the 
system. However, in 2D one is able to induce DW's in 
a magnetic layer up to a few thousand in a typical sam- 
ple, even for temperatures as low as a few K 8]. The 
use of this particular solution as a starting point is then 
justified, even at the low temperatures for which the ex- 
citations studied ahead are more important. 

The DW solution for Hh as obtained here bares resem- 
blance to the solitons studied in the ID Heisenberg chain 
0- As we see below, this similarity with the ID chain 
for the equilibrium solution does not lead to the same 
physics in 2D, the fluctuations about this solution are 
given by the usual massive magnons plus the additional 
gapless modes, only present in dimensions > 1. 

We consider small deviations from equilibrium given 
by 9 = 9(x) + £(r, t) and (f> = (f> + ri(r, t)/ sin9(x), where 
£,77 <C 1. Substitution of these relations in Eqs.© 
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where V(x) = (l/A 2 ) [2tanh 2 (x/X) - l] , and an identi- 
cal equation for rj(r, t). Dispersion relations can then be 
obtained from a standard choice of the solution in the 
form £(r, t) = b(r) exp — icut, which leads to the following 
eigenvalue equation, 



It is convenient to study the spectrum of the x depen- 
dent solutions to understand why there must be a gap 
in ID but not in 2D. This is understood if we realize 
that Eq.© yields the usual discrete eigenvalues associ- 
ated with the x direction and a mode that is given by just 
setting uj = 0. Let us call this zero eigenvalue solution 
by zero-mode for further reference. In ID, the zero-mode 
is a small static distortion from the original DW solu- 
tion localized near x = and with the functional form 
£(x) ~ sech(x/A), as it is seen by replacing x — *■ r in 
Eq.(|SJ) and integrating it for u> — 0. Since such a mode 
has no dynamics, the spectrum of the ferromagnetic ID 
anisotropic chain is gapped as it is well known. In par- 
ticular, away from the DW X 2 V(x) — ► 1, and a solution 
is promptly obtained both in ID and 2D, in the form 
of isotropic magnons with a massive spectrum given by 
w(k) = j(k 2 + l/A 2 ) where k is the wave vector. This 
dispersion holds also in 2D with k 2 = k 2 + k 2 . 

The key ingredient to distinguish 2D from ID is that 
with the additional degree of freedom, the zero-mode ac- 
quires dynamics and becomes the gapless magnon that 
we are studying. We see this by using a solution in the 
form b(r) = f(x)g(y) in Eq.JSJ), which separates in 



dx 2 



V(x) 



fix) 



/(*) 



and 



g"(y) = -Q 2 g(y)- 



(6) 



(7) 



[V 2 - V(x)] b(r) 



-b(r). 



(5) 



It is clear from these equations that the zero-mode associ- 
ated with the x direction in Eq.JfJJ is a propagating spin 
wave with dispersion given by U)(q) = jq 2 . In 2D this 
mode has a finite frequency and is localized about x — 0, 
as is seen from the solution b(r) = sech(x / X) exp iqy 
obtained from integrating Eq.© for ui(q) = jq 2 . An 
identical solution holds for rj(r,t). One can see that 
£ is related to transverse (_L y), while r\ to longitu- 
dinal oscillations. The complete solution is given by 
S(r,t) = S(r) + {dS(r)/d9)£(r,t) + yS r?(r,t), where 
£ and r\ are more generally linear combinations of the q 
and w dependent solutions worked out above. The gap- 
less dispersion corresponds to plane magnetization waves 
propagating with very small amplitudes inside the ID 
DW. 

It should be pointed that these are not the first gapless 
excitations reported in a 2D anisotropic chain. Gapless 
kinks that could propagate under very restricted condi- 
tions have been reported previously [lOj, and they con- 
sist basically of topological excitations inside the DW. 
Within our approximations, Eq.JJJ is linear and disper- 
sive and therefore will not sustain such localized kink-like 
solutions which will ultimately decay into ordinary spin 
waves. Gapless dispersions have also naturally appeared 
in solutions of the XXZ hamiltonian in 3D,|ll| showing 
that this is not a particularity of the 2D case we are 
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treating. Here it becomes clear the necessity of dimen- 
sions > 1 for the existence of such wall magnons. 

Consideration of these modes for the transport of elec- 
trons in magnetic materials appears to be absent from 
the up to date literature, where magnons are overwhelm- 
ingly considered in the analysis of measured data only at 
temperatures suffciciently high to be comparable to the 
anisotropy gap, whose typical values range from 50 K to 
150 K, depending on the material used. Absence of gap 
in the spectrum should make the contribution of these 
objects relevant for the resistivity at lower temperatures. 

In what follows, it will be seen that a transport theory 
for electrons interacting with the gapless magnons can in 
fact be put forward without complications and that the 
effects on the resistivity are expected to be quite observ- 
able for temperatures much lower than the anisotropy 
gap. 

Hydrodynamics is quantized by writing 



Z[qy] = 
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ensures that the con- 
in Eq.Q is written as 



and rj = £\qy — n/2]. This 
tribution of the free magnons 
7~(-m = Yl q ^qo\ a qi where a q is the (Heisenberg) bosonic 
operator that creates a magnon with wave vector q, 
huj q = JS 2 a 2 q 2 , and L y is the width of the sample. 

The electron- magnon interaction given by Eq.(J2J), can 
be put in the form Hi = Uq (x) + U\ (x) + H em > where 
Uq(x) = SJk&z tanh(x/ A), U\{x) — S Jxc 2: sech(x/A), 
and 
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The term given by U(x) = Uq{x) + U\(x) affects the 
electron's motion only in the x direction and can be in- 
corporated as an effective potential in a purelly electronic 
hamiltonian, H e — — h 2 V 2 /2m+U(x). With these defini- 
tions, the full hamiltonian IC + Hh + Hi can be rewritten 
as H e + H m + Hem j where effectively, the two first oper- 
ators are single-particle hamiltonians while the last one 
accounts for the interaction. 

As a result of its interaction with the DW, the electron 
experiences an effective potential that combines a spin- 
conserving term, Uq{x) oc a z and a spin- flipping term 
U\ (x) oc a x . The former has the approximate form of a 
smooth step barrier centered at the DW while the latter 
exists only within the DW. From this structure of U(x) 
it is seen that spin-up electrons tend to be blocked from 
left to right while the opposite is true for the spin-down 
electrons. This is expected since the energy of one spin 
rises when it reaches the region where the magnetization 
is opposed to it. The DW tends to work then as a spin 
diode, since it will conduct one species more favorably in 
one direction. 

Interaction with magnons, given by H e m, may be 
viewed as a complex potential in the x direction, that 



besides flipping the spin of the electrons, has also a dis- 
sipative character. Spin-flip makes transport favorable 
since the electron can travel into the opposite-spin re- 
gion with an energy gain. However, the momentum lost 
by emitting a magnon causes the electrons to loose energy 
in the process, and this is a temperature dependent com- 
ponent. What we expect, on a qualitative basis, is that 
these magnons will lower the efficiency of such a diode 
as a spin filter, having the same effect on the transport 
electrons as the usual hot magnons, but showing their 
effect at much lower temperatures. 

To calculate the conductivity due to such an interac- 
tion we start with the Boltzmann transport equation, 
which may be written in a relaxation time approximation 
by assuming that the change in the electrons' distribution 
function is given by (5/ CT (k) = C(fc)k • E, where E is an 
in-plane uniform electric field perpendicular to the DW. 
To leading order, the collisions with the DW affect only 
the scatering ampliudes W£ k , = W{£, + W^, a , where 
Wkk' ^ S ^ ne probability per unit time for an electron in 
an initial state |ktr) to be scattered to |k'(r'). With these 
assumptions, the relaxation time obtained is 
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where A is the area of the sample and 9 is the angle 
between k and E. The crossed term that appears on the 
right side of Eq.JSJ) is a particularity of the 2D system, 
it is absent in both 1 and 3 dimensions. It does not 
modify the relaxation time considerably though since, as 
it is known, the collisions that mostly affect the transport 
current consist of nearly backward processes (k ~ — k') 
with exchange of momentum ~ 2fci?. The conductivity is 
readily calculated from r cr (k), 
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The scattering amplitudes are calculated by using the 
exact states of K. + Uq(x) + 7Y m , which is the only 
part in H that does not vanish far from the DW. 
If we denote such states by |kcrng), where n q indi- 
cates the magnon state in occupation-number formal- 
ism, then the amplitudes will be given by W££, = 
2TrN/hJ2 qq > KWqlV^'lkVrv)! 2 ^, where N is the 
number of domain walls and Sq constrains the transitions 
within the ones that conserve total energy and momen- 
tum. Here, V 7 " 7 = V X + V y and V aa = V x , while Vi is 
the sum of all terms in U\ (x) + H em that are oc Ui . As 
stated, the detailed solutions constitute text book mate- 
rial, but are nonetheless sufficiently awkward not to be 
shown here^2- The solutions are used to numerically 
handle Eqs.JSJ and © in obtaining the resistivity p. 

Fig [21 compares the resistivity along the x direction of 
a 2D sample for different values of the temperature as 
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FIG. 2: Resistivity in the x direction as a function of fcfA 
where A is the size of the DW. The continuous line corresponds 
to the elastic case (no magnons) , while the dotted and dashed 
lines were computed including the electron-magnon interac- 
tion for T=10 K and 25 K respectively. The anisotropy gap 
used to obtain these curves was about 116 K (10 meV) and 
in all cases 3k = 0.5 eV. The inset shows the temperature 
dependence for /cfA ~ 4.2. 

a function of kp\ with and without including the gap- 
less magnons. The sample size has been taken as L x = 
500 /im and L y = 50 /im, and N ~ 3000. Such a large 
number of DW are not expected to be naturally pesent 
in magnetic layers down to arbitrarily low temperatures. 
However, the interesting nanolayers for spintronics and 
device applications should present DW at much lower 
temperatures, and such a number is achieved by exper- 
imental techniques down to 2 K[8j. If one considers the 
typical electronic densities, in the range from 10 15 to 10 20 
cm -2 , in layers, the DW sizes shown in Fig|2] range from 
a few A to several tenths of A, consistent with what is 
also observed in experiments 8] . 

One sees that an observable increase in the resistivity 
results for most values of fcpA, as expected due to the 
localized character of the gapless magnons. What is no- 
table here, is that even at a temperature of 10 K the effect 
can be quite observable, while the usual gapped magnons 
considered in most of the literature would, in this case, 
affect the results only for temperatures larger than about 
116 K (the anisotropy gap). As a final remark, we would 
like to point out that these effects can be experimentally 
checked by measuring the magnetoresitance related to 
fields higher than the critical DW field. This eliminates 
scattering processes as phonons or non-magnetic impuri- 
ties from the background, as it is known. If the remain- 



ing resistivity shows a temperature dependence consis- 
tent with a magnon bath for temperatures considerably 
lower than the anisotropy gap, this will strongly suggest 
the presence of the gapless magnons. 

In closing, we have found that, contrary to the case 
in ID and to common belief, the magnetic excitation 
spectrum of a 2D (single-layer) anisotropic ferromagnet 
should influence charge transport at considerably low 
temperatures in spite of the anisotropy. This is due to 
the existence of one dimensional gapless magnons that 
can propagate in the domain walls formed as the lowest 
energy configuration. We have provided a theory to de- 
scribe how electrons interact with both these magnons 
and the domain wall. We have presented results showing 
the effect of these magnons on the transport of electrons 
through domain walls. 
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